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ABSTRACT 

We investigate with three-dimensional numerical simulations of supernova 
driven turbulence in the interstellar medium (ISM) the relationship between the 
velocity dispersion of the gas and the supernova rate and feedback efficiency. 
Our simulations aim to explore the constancy of the velocity dispersion profiles 
in the outer parts of galactic disk at ~ 6 — 8 km s _1 , and the transition to 
the starburst regime i.e., high star formation rates associated with high velocity 
dispersions. With our fiducial value of the supernova feedback efficiency (i.e., 
e = 0.25 corresponding to an injected energy per supernova of 0.25 x 10 51 ergs), 
our results show that (a) supernova driving leads to constant velocity dispersions 
of a ~ 6 km s" 1 for the total gas and am ~ 3 km s" 1 for the Hi gas, independent 
of the supernova rate, for values of the rate between 0.01 and 0.5 the Galactic 
value (t]g), (b) the position of the transition to the starburst regime (i.e., location 
of sharp increase in the velocity dispersion) at around SFR/Area ~ 5 x 10~ 3 — 
10~ 2 M yr _1 kpc~ 2 observed in the simulations, is in good agreement with the 
transition to the starburst regime in the observations (e.g., NGC 628 and NGC 
6949), (c) for the high SN rates, no Hi gas is present in the simulations box, 
however, for the total gas velocity dispersion, there is good agreement between 
the models and the observations, (d) at the intermediate SN rates (rj/r] G ~ 
0.5 — 1), taking into account the thermal broadening of the Hi line helps reach a 
good agreement in that regime between the models and the observations, (e) for 
v/Vg < 0.5, a and cxm fall below the observed values by a factor of ~ 2. However, 
a set of simulation with different values of e indicates that for larger values of the 
supernova feedback efficiencies, velocity dispersions of the Hi gas of the order 
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of 5 — 6 km s -1 can be obtained, in closer agreement with the observations. 
The fact that for r]/r] G < 0.5, the Hi gas velocity dispersions are a factor ~ 2 
smaller than the observed values could result from the fact that we might have 
underestimated the supernova feedback efficiency. On the other hand, it might 
also be an indication that other physical processes couple to the stellar feedback 
in order to produce the observed level of turbulence in galactic disks. 

Subject headings: Galaxies — Galaxies : velocity dispersion-ISM : supernova- 
turbulence-instabilities: thermal instability 

1. INTRODUCTION 

It is now a well established fact that the interstellar medium (ISM) in galactic disks 
is turbulent (Larson 1981; Scalo 1987; Dickey & Lockman 1990; Elmegreen & Scalo 2004). 
In most spiral galaxies, and after correcting for instrumental effects, the vertical velocity 
dispersion derived essentially from H I observations is observed to vary radially from ~ 12—15 
km s _1 in the central parts to ~ 4 — 6 km s~ x in the outer parts (van de Kruit & Shostak 
1982 for NGC 3938; Shostak & van der Kruit 1984 for NGC 628; Dickey, Hanson & Helou 
1990 for NGC 1058; Kamphuis & Sancisi 1993 for NGC 6964; Rownd, Dickey & Helou 
1994, for NGC 5474; Meurer, Mackie & Carignan 1994 and Meurer et al. 1996 for NGC 
2915; de Blok & Walter 2005 for NGC 6822). In most cases, the latter values exceed those 
expected by the thermal broadening of the Hi emission line. As an example, the radial 
dependence of the velocity dispersion in NGC 1058 is shown in Fig. 1. Fig. 2 shows the 
behavior of the characteristic velocity dispersion of a sample of galaxies when plotted versus 
their characteristic star formation rate. The galaxies shown in Fig. 2 correspond to a sample 
where both the velocity dispersion and the star formation rate were available in the literature. 
Similarly to the observed velocity dispersion radial profiles, the data plotted in Fig. 2 argues 
for the existence of a minimum level of turbulence in galactic disks. The case of NGC 2915 
is particularly intriguing as it has an extended H I disk with a constant velocity dispersion, 
even in the outer parts, where stars do not form (Meurer et al. 1996). 

Many physical processes might contribute to the driving of the turbulence in the ISM, 
acting on different scales and injecting different amounts of kinetic energy into the medium. 
These driving mechanism can be related to stellar activity (e.g., ionizing radiation from 
Hn regions, jets from young stellar objects, stellar winds, supernova explosions), or to gas 
hydrodynamical or magnetohydrodynamical instabilities (e.g., thermal and gravitational in- 
stabilities, magnetorotational instability, Kelvin-Helmotz, Rayleigh- Taylor and Parker in- 
stabilities) or, eventually, to the complex interaction between these processes. Mac Low & 
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Klessen (2004) showed, using simple analytical estimates, that the global energy input into 
the ISM of the Galaxy from supernova explosions (SNe) is much larger than from any of the 
above mentioned processes. However, their estimates are global ones which do not compare 
the efficiency of each process at different galactic radii and the possible interactions between 
two or several processes. Other, yet poorly explored, sources of kinetic energy injection into 
the ISM of some galaxies might reside in perturbations occurring on galactic scales (e.g., 
frequent minor mergers, tidal interactions with satellite galaxies and ram pressure effects). 

In the past decade, some of the above mentioned energy injection mechanisms into 
the ISM have been investigated by means of numerical simulations. Wada & Norman 
(1999,2001), Wada, Spaans & Kim (2000) and Wada, Meurer & Norman (2002) investigated 
the evolution of thermal and gravitational instabilities (TI and GI, respectively) in galactic 
disks. Wada, Meurer & Norman (2002) showed, using high resolution two-dimensional sim- 
ulations of an NGC 2915 like disk that the turbulent energy spectra can be maintained in a 
quasi- stationary state in the absence of stellar feedback. In the latter simulations, the kinetic 
energy decay and radiative cooling are compensated by the interaction of the galactic shear 
with the gas self-gravity. Another indication that energy might be injected into the ISM of 
low star forming galaxies on large scales before cascading to smaller scales via thermal and 
gravitational instabilities is provided by the analysis of their H I gas morphology. By compar- 
ing the Hi emission map of Holmberg II (Ho II) to synthetic Hi maps of driven turbulence 
which include also cooling, heating and gravity, Dib & Burkert (2005) found that kinetic 
energy is injected into the ISM of Ho II on a scale of ~ 6 kpc, which is much larger than 
the scale implied by SN driving. Galactic shear might play a role in supporting turbulence 
in Ho II, however, the comet shape structure of the Hi gas in Ho II (Bureau et al. 2004) 
suggests that turbulence might be also partially induced by the effect of tidal interactions 
with Ho II's two satellite galaxies. 

The coupling of the galactic shear to magnetic fields can trigger the magnetorotational 
instability (MRI) (Balbus & Hawley 1991). MRI has been the scenario invoked by Selwood 
& Balbus (1999) to explain the constancy of the velocity dispersion in the outer parts of 
NGC 1058. The main argument of Selwood & Balbus (1999) against a SN-driven ISM 
in the outer regions of NGC 1058, in addition to the fact that column densities are low 
and star formation is inefficient, is the observed uniformity of the Hi velocity dispersions 
(Dickey, Hanson & Helou 1990) whereas Hn regions in NGC 1058 are primarily observed in 
narrow spiral arms (Ferguson et al. 1998). Local, three-dimensional simulations with a single 
phase medium by Kim, Ostriker & Stone (2003) of the MRI result in velocity dispersions 
of the order ~ 1.6 — 3.2 km s _1 . Global one-phase medium simulations by Dziourkevitch 
et al. (2004) and Dziourkevitch (2005) show that some amount of turbulent motions can 
be created by the MRI in galactic disks. However, this turbulence is mostly located in the 
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inner parts of the disk (the inner 4 kpc) and the velocity dispersion does not exceed in that 
case ~ 3.5 — 4 km s _1 , dropping quickly to very small values at larger radii. Two and three- 
dimensional MRI simulations in a medium affected by the TI but with no stellar feedback 
have been presented by Piontek & Ostriker (2004,2005). In the three-dimensional models, 
the latter authors show that the velocity dispersion depends on the average density of the 
medium (larger dispersions at lower densities). At the lowest density they have considered 
(0.25 cm -3 ), the one-dimensional velocity dispersion is found to be 2.2, 2.2 and 1.1 km s -1 
for the three directions of the box, respectively (~ 4.5 km s" 1 for the three components if 
thermal broadening is not subtracted). The latter values are a factor of 3-6 lower than the 
observational values (~ 1.5-2 if thermal broadening is not subtracted). Recent simulations of 
the MRI in a multiphase medium with star formation included show that the MRI might be 
completely suppressed/overwhelmed by stellar feedback (M. Korpi, private communication). 

Another major source of energy input into the ISM is the stellar energy feedback, par- 
ticularly from massive stars. The latter can be delivered to the ISM in the form of ionizing 
radiation and stellar winds from O and B stars (Kessel-Deynet & Burkert 2003) and from 
clustered or field supernova explosions (McKee & Ostriker 1977). Two and three-dimensional 
numerical models of SN explosions models in the ISM have been presented in the literature 
(Rosen & Bregman 1995; Korpi et al. 1999a,1999b; Gazol-Patiho & Passot 1999; de Avillez 
2000; de Avillez & Berry 2001; Avila-Reese & Vazquez-Semadeni 2001; Kim et al. 2001; 
Kim 2004; de Avillez & Breitschwerdt 2004,2005; Slyz et al. 2005; Mac Low et al. 2005). 
In these models, the authors have focused on problems like the evolution of SNe bubbles 
and their outburst through the galactic disk, the halo-disc interaction, the vertical scale 
heights and volume filling factors of the different gaseous phases and the effects of super- 
nova explosions on the Galactic dynamo (Ferriere 1992a,1992b, 1998a, 1998b, 2000). However, 
the existence of a correlation between the SN rate and the velocity dispersion of the gas 
has not been investigated so far. A different interpretation of the velocity dispersion in the 
outer regions of galactic disks is discussed by Schaye (2004), which argues that if the sur- 
face density in the outer parts of the disk falls below a critical value (i.e., column density 
Log(N H ) < 27.5 cm -2 ), the medium can be efficiently heated to a temperature of T ~ 8000 
K by a background of ultraviolet (UV) ionizing radiation, both of galactic and extragalactic 
origin. Thus, thermal broadening could account for the observed level of line broadening in 
the outer parts of galactic disks, leaving little room for a dynamical, turbulent component. 
This hypothesis can be appropriately tested if the level of UV radiation, particularly the 
extragalactic component, and the temperature radial profiles in galaxies would be better 
constrained from observations. 

The aim of this paper is to assess how much of the velocity dispersion observed in the 
ISM of galaxies can be due to SN feedback for various values of the SN rate and feedback 
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efficiency. A particular point of interest is to understand the constancy of the velocity 
dispersion in galactic disks at different radii where the star formation rate is expected to 
decrease with increasing radius from the galactic center as predicted by many empirical 
star formation laws (Schmidt 1959,1963; Kennicutt 1998a, 1998b; Dopita & Ryder 1994; 
Prantzos & Silk 1998; see also Li et al. 2005a,2005b who finds, as an explanation of the 
global Schmidt laws observed in galaxies, that the star formation rates correlate with the 
efficiency of gravitational instability). This paper is organized as follows. In § 2, we describe 
our models and the relevant parameters. In § 3, we describe how synthetic observations 
are derived. The velocity dispersion dependence on the feedback efficiency, supernova rate 
and average gas number density is presented and discussed in § 4, § 5 and § 6 respectively. 
Detailed comparison to the observations is performed in § 7. In § 8 the need for improved 
numerical models is critically reviewed and in § 9, we summarize our results and conclude. 



2. THE MODEL 

In order to understand how the SN rate and energy feedback efficiency affect the ve- 
locity dispersion of the gas, we resort to a simple numerical model in which the vertical 
stratification, galactic rotation, magnetic fields and the gas self-gravity are not included. We 
use the ZEUS-3D code (Stone & Norman 1992a,b) to solve the equations (mass, momentum 
and internal energy conservation) of ideal gas dynamics. We simulate a 1 kpc 3 volume of 
the ISM with a grid resolution of 128 3 . Periodic boundary conditions are imposed in the 
three directions. In most simulations and if not specified otherwise, the initial density field 
is homogeneous with a number density of n = 0.5 cm~ 3 . The velocity and temperature are 
everywhere zero and 10 4 K, respectively. Radiative cooling of the gas is included by directly 
interpolating in the the solar metallicity cooling curves of Dalgarno & McRay (1972) in the 
temperature range of [100 K,10 4 K] and of the more recent data of Sutherland & Dopita 
(1993) for the temperature range of [10 4 K-10 8,5 K]. The gas is not allowed to cool below 
the minimum temperature T min = 100 K. The maximum temperatures at the explosion sites 
reach values of ~ 60 — 70 x 10 6 K which makes the hot gas fall in the hot stable regime. 
A polytropic equation of state with a specific heat ratio of 5/3 is used. This proves to be 
justified because of the lower cutoff temperature of 100 K that limits over-densities to a few 
tens of cm~ 3 and therefore the gas remains mostly monoatomic. 

Stellar feedback is modeled as resulting from type II SN explosions only and the energy 
is injected into the ISM instantaneously. The total energy of each explosion is taken to be 
E SN = 10 51 erg (Chevalier 1977; Abbott 1982; Woosley & Weaver 1986 and Heiles 1987). 
Only a fraction of the total energy is transferred to the ISM in the form of thermal energy. 
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This defines the feedback efficiency parameter e. The energy is injected into 3 3 cells around 
the central explosion cell following a Gaussian profile. The energy each neighboring cell 
receives is weighted by the square of its distance to the central cell in order to insure a better 
isotropy of the explosion. The site where a new SN explosion occurs is chosen randomly 
under the condition that the local density n is such that n > n. This assumption leads to a 
more realistic fraction of clustered SN explosions. The time interval between two consecutive 
SN explosions is given by At$N = 1/^7, where 77 is the SN explosion rate. However, when 
the time interval between two consecutive SN explosions become shorter than the CFL 
(Courant-Friedrisch-Levy) time step, more than one SNe are detonated simultaneously at 
different locations of the grid. The number of SNe detonated in that case is taken to be the 
closest integer to the ratio dt/AtsN, where dt is the CFL time step. We define a Galactic 
SN rate (rjc) of 2.58 xlO -4 yr _1 kpc~ 3 assuming a Galactic radius of 15.5 kpc and a scale 
height for type II SN of 90 pc (Miller & Scalo 1979). We use a frequency of 1/57 yr -1 
(Capellaro, Evans & Turatto 1999) which is smaller than earlier estimates of 1/50 yr _1 by 
van den Bergh & McClure (1990) and 1/37 yr _1 by Tammann, Loffler & Schroder (1994). 
Note that in our simulations, heating by a UV background radiation is not accounted for. 
It has not been included because no clear recipe exists that correlates the amount of UV 
heating that should be added with the different SN rates. One possibility is to distribute a 
certain fraction of the energy associated with a given SN rate on the total grid. However, at 
low SN rates, associated with low density environments in the outer parts of galactic disks, 
the disk might be more easily heated by incident UV photons of extragalactic origin which 
intensity is not well determined (see discussion in Schaye 2004). 

In order to test the effects of numerical resolution, we explode a single supernova in a 
medium of initial temperature of 10 4 K, initial homogeneous and uniform density of 0.5 cm~ 3 
on a grid representing a physical scale of 1 kpc. The test is similar to the one presented by 
Mac Low et al. (2005). However the latter authors, aside from having a smaller simulations 
box (200 pc), assumed the SN remnant to evolve in a medium of negligible pressure (i.e., 
temperature of 10 K and average density of 0.1 cm -3 ). Fig. 3 displays the time evolution of 
the shell position R s h, measured as being the position of the density peak, in simulations with 
numerical resolutions of 64 3 , the fiducial resolution of 128 3 , and 256 3 . Clearly, a resolution 
of 64 3 is not sufficient to resolve properly the dynamics of a single SN remnant. On the 
other hand, particularly after the first ~ 0.1 Myrs in the lifetime of the SN remnant, the 
discrepancy in the position of the shell at the resolutions of 128 3 and 256 3 is of the order 
of 10 % and does not exceed the value of 20 %. We conclude that our fiducial resolution of 
128 3 is good enough for the purpose of studying the global dynamical effects of the energetic 
input of SN explosions into the ISM. 
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3. 



ANALYSIS AND DERIVATION OF THE OBSERVABLES 



The simulations are evolved until kinetic energy reaches a stationary value. Fig. 4 shows 
the evolution of kinetic energy in a number of simulations. The equilibrium value for the 
kinetic energy is reached when the dissipation equals the amount of injected kinetic energy. 
The medium acquires kinetic energy from the SNe explosion-induced thermal pressure gra- 
dients, and the thermal pressure gradients associated with TI which occurs in the dense 
expanding shells. In most simulations, the equilibrium of thermal energy is also reached 
except for the simulations with a high supernova rate (> 2.5 r/c)- In the latter simulations, 
the overlap radii between supernova remnants is very small, and all the gas continues to heat 
up as more and more energy is injected into the system. Fig. 5 and Fig. 6 show snapshot 
two-dimensional cuts for models with (i]/i]G,€) = (0.1,0.25) and (1,0.25), respectively. In 
the simulations with the lower rates, larger and denser clouds are able to form under the 
effect of TI, before being eventually dispersed by a next generation local SN explosion. 

We calculate the velocity dispersion of the gas in two complementary ways. In the first 
method, we evaluate the characteristic mass weighted velocity v c from the three-dimensional 
data as 



where the index i runs over the number of cells in the simulation box. An average 
is made over the last 5 Myrs (5 values) in each simulations in order to smooth for time 
fluctuations which is particularly useful for the low SN rate simulations. On the other 
hand, following a more observational approach, the one-dimensional velocity dispersion a is 
obtained by fitting the mass-weighted line of sight velocity profile. Intensity is assumed to 
be proportional to the mass along the line of sight. This is particularly true in the case of 
the Hi line (Rholfs & Wilson 1996). The velocity profile is then normalized to its maximum 
value. A second velocity dispersion which we call am is obtained by fitting a velocity profile 
where only cells which have temperatures < 12000 K and number densities n > 0.25 cm -3 
have been accounted for, thus mimicking the velocity profile of an H I emission line. We have 
tested the dependence of the fit parameters on the size of the adopted velocity bin. The 
relevant parameter (i.e., width of the velocity profile) is practically unchanged as we vary 
the velocity bin size from 0.1 km s -1 up to 2.5 km s -1 , only the fit-error on the parameters 
changes, remaining however very close to the value of the spectral resolution. The results we 
will show correspond to a spectral resolution of 1 km s _1 . The latter value is characteristic 
of single-dish radio telescopes (e.g., Effelsberg radio telescope, Green Bank radio telescope) 
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and nears the spectral resolutions obtained with the VLA (Very Large Array) ~ 2.5 km s -1 
which will be enhanced when the EVLA (Extended Very Large Array) becomes operational. 
Here also, for each simulation, an averaging over the last 5 Myrs has been performed (5 
estimates). Errors on a and o"hi are average values of the individual errors derived from the 
parameters of the fit functions whereas the error on v c is simply a statistical error over the 
5 estimates. 

We have attempted to fit the line of sight velocity profiles with several functional distri- 
butions, namely, a Gaussian, a Lorentzian, a Moffat profile (modified Gaussian) and a Voigt 
profile (see Lang 1980 for the mathematical definition of each function). In all cases, the 
Gaussian, Lorentzian and Moffat fit functions yield the same dispersion value which has a 
slight different meaning in each case (see Lang 1980 for details). Fitting with the Voigt pro- 
file works only in a limited number of cases, but in those case proves to be a better fit of the 
profiles. Examples are shown in Figs. 7-10. Fig. 7 and Fig. 8 show the 1 km s" 1 binned total 
gas and Hi gas line profiles for the model with (t]/t]g,£) = (1,0.25), respectively, whereas 
Fig. 9 and Fig. 10 show the same profiles for the model with {r//r]G, e) = (0.1, 0.25). Though 
Gaussian fitting proves to be quite satisfactory, it is worth mentioning at this stage that 
our simulated velocity profile have wings that are slightly broader than those associated to 
Gaussian functions. This non-Gaussianity of the line profiles has been already pointed out 
by Ricotti & Ferrara (2002) in their Monte-Carlo models of the ISM dynamics. In a realistic 
and turbulent ISM, one expects that the turbulent contribution to the line broadening would 
be described with a Lorentzian and the thermal broadening of the line by a Gaussian. The 
convolution of the Gaussian and Lorentzian function results in a Voigt profile. However, 
even when we introduce a thermal broadening component to the line profile (in § 7), the 
Voigt profile fails to systematically fit the line profiles of the different models. 

Since energy injection in our simulations is discrete, a similar value r\ x e with different 
permutations of rj and e does not yield similar results because of the differences in the inter- 
actions of the expanding shells and non-linear development of TI. Ideally, a full investigation 
of the two-dimensional parameter space would be necessary, however, this would lead us 
beyond our current computational capabilities. At this stage, we have taken an intermediate 
approach and have explored the effects of i] and e independently by fixing one parameter 
and varying the other. 



4. THE EFFECT OF THE FEEDBACK EFFICIENCY 



We performed a first set of simulations for which we varied the supernova feedback 
efficiency between 0.05 and 1 for a constant value of the supernova rate, 77 = 0. 1 r\ G . Fig. 11 
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shows the dependence of the characteristic velocity v c , line of sight velocity dispersion a 
and Hi line of sight velocity dispersion <7hi on e. Assuming that all of the three values 
would be zero if the feedback efficiency is zero, this would mean that v c , a and cr H i would 
rapidly increase with increasing e, followed thereafter by a slower increase at larger values 
of e. As a velocity has the dimensions of the square root of an energy, we have attempted 
to fit v c , a and <7hi by functions of the form Ay/e. Fits with these functional forms are 
over-plotted to the data in Fig. 11. The values of the fit parameter A are 17.15 ± 0.59 km 
s -1 , 7.91 ± 0.36 km s" 1 and 6.49 ± 0.44 km s -1 for v c , a and <7 H i, respectively. In all three 
curves the simulations are noticeably higher than the estimate of the fit functions for the 
small values of e. This is plausibly an indication that a second mechanism is responsible for 
generating kinetic energy in the medium on top of the direct supernova driving in the regime 
of low energy injection. TI is likely to play an important role in that regime as clouds have 
time to condense further, thus enhancing the effects of TI, before being destroyed by a local, 
next generation SN explosion. We will particularly focus on the role of TI in the low energy 
injection regime in the next section in which we investigate the dependence of the velocity 
dispersion on the SN rate and for which we have more available simulations. 

The total feedback energy (kinetic+thermal) deposited in the ISM by a type II SN 
exploding in a medium of average density 1 cm" 3 and solar metallicity was estimated by 
Thornton et al.(1998) to be ~ 7 x 10 50 ergs at the time of the peak luminosity (i.e., t ), 
dropping to ~ 0.2 x 10 50 ergs at t ~ 13i , roughly equally distributed in thermal and 
kinetic energy. In our models, t coincides with the time at which the energy is released 
for each SN explosion and distributed in a volume of radius ~ 11.7 pc. However, the 
estimates of Thornton et al. (1998) should be regarded with some caution essentially because 
they are based on a one-dimensional model in which the gas has few degrees of freedom. 
Furthermore, the cooling curve adopted by Thornton et al. (1998) excludes the radiative 
cooling from neutral atoms for temperatures lower than 10 4 K and does not include other 
potential cooling mechanisms such as neutrino cooling, bremsstrahlung radiation and cooling 
by thermal conduction. For the investigation of the supernova rate effect on the velocity 
dispersion, we adopt a conservative value of the feedback efficiency of e = 0.25. 



5. THE SUPERNOVA RATE- VELOCITY DISPERSION RELATION 

The most relevant parameter in our simulations is the supernova rate rj. Fixing e at a 
value of 0.25, we performed a set of simulations with different values of r] ranging from 0.01 
to 10 times the Galactic value. The dependence of v c and a on r\ is displayed in Fig. 12 
and Fig. 13, respectively. For rj/rjo > 0.5, Fig. 12 and Fig. 13 show that v c and a increase 
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rapidly with an increasing supernova rate, thus mimicking a starburst regime similar to the 
one observed in Fig. 2. A fit to the data for r]/r)o between 1 and 10 is a = 0.78(±0.26) i]/r)G + 
8.30(±1.36) km s _1 . This relation might be useful for semi-analytical modeling of the central 
regions in galaxies, e.g., the modeling of active galactic nuclei (AGNs), in which hydrostatic 
equilibrium requires that the effective gravitational potential be balanced by the turbulent 
pressure of the gas. In current analytical models of AGNs, the intrinsic turbulent velocity 
of the gas is neglected and gas clouds are assumed to have the same velocity dispersion 
of the nuclear stellar cluster (e.g., Schartmann et al. 2005). For r]/r] G < 0.5, the velocity 
dispersion shows a slow decrease with a decreasing rj. This transition at r]/r] G — 0.5 could 
indicate that the velocity dispersion of the medium in the low rate regime is not fixed by SN 
driving alone. SNe explosions will cause a certain fraction of the gas to be maintained in the 
thermally unstable regime when cold gas is restored to the warm phase. Thermal instability 
occurs in the cooling expanding shells, but also everywhere in the inter-supernova remnants 
medium where the criterion for thermal instability (see Eq. 9 in Vazquez-Semadeni, Gazol 
& Scalo 2000) is satisfied, thus adding an extra component to the kinetic energy injected 
into the medium. Since we are using the same cooling curve in the low temperature regime 
as Vazquez-Semadeni et al. (2000), the fit used by the latter authors and by Spitzer (1978) 
for the data of Dalgarno & McRay (1972) remains valid and the thermally unstable regime 
will be confined in the temperature range 398 K < T < 10000 K. 

We interpret the flatness of the a — rj relation for r]/r] G < 0.5 as resulting from the 
interplay between direct supernova driving and thermal instability. Fig. 14 shows, after con- 
vergence is reached, the dependence of the volume filling factor of the thermally-unstable gas 
on the normalized SN rate, whereas Fig. 15 and Fig. 16 display, as examples, the time evolu- 
tion of the volume filling factor of the unstable (398 K < T < 10000 K), cold (T < 398 K) and 
warm gas (T > 10000 K) in the models with {r]/r] G , e) = (0.01, 0.25) and (77/770, e) = (0.1, 0.25), 
respectively. In Fig. 14, the volume filling factor of the unstable gas which is an indicator 
of the occurence of TI shows a transition at i]/i] G ~ 0.5 which roughly corresponds to the 
position of the transition observed in the r) — v c and r] — a relations, and a non-zero value 
at smaller values of rj. For the large SN rate values, the gas is predominantly hot, becoming 
increasingly hotter with time. Therefore the volume filling factor of the unstable gas in that 
regime is close to zero. In the regime where rj/rjo < 0.5, TI is more efficient in converting 
the gas into the cold phase at the lower end values (i.e., TI is more efficient at i]/i] G = 0.01 
than at i]/i] G = 0.1). Fig. 15 and Fig. 16 show that the converged value of the volume filling 
factor of the cold gas in the simulation with (rj/rja, e) = (0.01, 0.25) is F c ~ 0.8, whereas this 
value is F c ~ 0.5 for the simulation with (rj/rja, e) = (0.1, 0.25). For the same average density 
in both simulations, SNe exploding in a medium with larger fractions of its volume in the 
cold phase as in the simulation with {r\jr\ Gl e) = (0.01, 0.25) will evolve in a cooled lower pres- 
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sure environment than in the simulation with (rj/rja, e) = (0.1, 0.25). Numerical simulations 
of SN driven turbulence by Kim (2004) show that the evolution of SNe remnants in media 
with lower external pressures leads to higher velocity dispersions of the gas. The existence 
of a background heating process could modify to some extent the present conclusion. How- 
ever, the background heating should be strong enough to oppose the dramatic cooling of 
the medium in the low SN regime and help maintain large fractions of the gas at higher 
temperatures. However, in the absence of an external heating mechanism to the galaxy, the 
backgournd heating can only be a small fraction of the energy released by SN explosions and 
might not play a significant role (see also the discussion in § 8). Finally, it is clear that for 
very small SN rates, SN driving would not be able to sustain any turbulence in the medium. 
Turbulence will then decay before the next SN explosion occurs. 



6. THE EFFECT OF THE AVERAGE DENSITY 

As the simulations described in this chapter are not scale free, because of the presence 
of a realistic cooling function, another relevant point to investigate is the role of the average 
density. Empirical star formation laws state that the star formation rate decreases with 
decreasing gas surface density (Schmidt 1959,1963; Kennicutt 1998a, 1998b; Dopita & Ryder 
1994; Prantzos & Silk 1998). In all previous simulations, we have used an average density 
of 0.5 cm -3 . It is particularly interesting to test the effect of varying the density for small 
values of the supernova rate. Over-plotted on Fig. 13 are the results of two simulations with 
v/Vg = 0.05 and rj/r/c = 0.01 where the average density have been decreased, by a factor of 
5 and 10, respectively. These are shown with the full square and full hexagon, respectively. 
Dropping the scaling coefficient, if one assumes the supernova rate (i.e., star formation rate)- 
gas density to follow a Kennicutt type law r]/rj G = n 1A±0 - 15 (Kennicutt 1998a, 1998b), for an 
average density value n = 0.5 cm~ 3 corresponds rj/r/c = 0.38±[j;Q37. To values of rj/r/G = 0.05 
and 0.01 will correspond, using the same law, average densities of n = 0.117±oq26 cm -3 and 
n = 0.037±oqJ23 cm -3 , which roughly equal the densities of 0.1 cm -1 and 0.05 cm -1 we have 
adopted for those rates. Thus, we can consider the three points (77/77^, n[cm~ 3 ]) = (0.5, 0.5), 
(0.05,0.1) and (0.01,0.05) as being a rough representation of a Kennicutt law in the 77 — a 
space. Tentatively, the preliminary conclusion we can draw here is that the flatness of the 
a — r] relation around ~ 6 km s" 1 can be maintained at low values of the supernova rate if 
the average density is reduced for lower rates, as predicted by the Kennicutt law. However, 
more simulations are needed to confirm this result and to probe the results for other star 
formation laws. Unfortunately, this is beyond the scope of this work, essentially for reasons 
of CPU time. 



7. COMPARISON TO THE OBSERVATIONS 



In Fig. 17, the same simulations appearing in Fig. 13 are shown after transforming the 
SN rate per unit volume into a star formation rate per unit area (units of M yr _1 kpc~ 2 ). 
We use the transformation of the SN rate into a star formation rate (SFR), r]/SFR=0.0067, 
derived using the PEGASE stellar population synthesis model (Fioc & Rocca-Volmerange 
1997) and assuming a Salpeter Initial Mass Function (IMF) (Salpeter 1955). We perform a 
comparison to two galaxies, NGC 628 and NGC 6946 for which the star formation rates have 
been estimated, at different radii, from Ha observations (Martin & Kennicutt 2001), along 
with velocity dispersion estimates which are derived from H I 21 cm line observations (Shostak 
& van der Kruit 1984 and Kamphuis & Sancisi 1993). Fig. 18 shows a comparison of the 
velocity dispersion measured from the Hi gas (100 K < T < 12000 K) velocity profile <thi to 
the same observations. A number of remarks can be drawn from the comparisons presented in 
Fig. 17 and Fig. 18 : (a) The position of the transition to the starburst regime (i.e., location 
of sharp increase in the velocity dispersion) at around SFR/Area ~5x 10~ 3 — 10~ 2 M Q yr _1 
kpc -2 observed in the simulations, is relatively in good agreement with the observations 
in NGC 628 and NGC 6949. It is also in very good agreement with the transition to the 
starburst regime observed in Fig. 2, (b) there is good agreement between our models and the 
observations at the high SN rate values within 2 — 3 km s -1 . This difference can be easily 
explained by the effect of beam smearing which tends to increase the observed velocity 
dispersions, particularly in the inner parts of galaxies (c) At intermediate and low SN rates, 
a and <7hi fall below the observed values by a factor of 2-3 even when the density correction 
related to the Kennicutt law is taken into account in the case of crm (open square and star 
in Fig. 18). The values of a and <t H i become very similar at the low rates as most of the 
gas in these simulations has temperatures that are below 12000 K (see the profiles in Fig. 9 
and Fig. 10). On the other hand no Hi gas is found in the simulations with rj/rjo > 1- It 
should however be kept in mind that we have adopted a rather conservative value for the 
supernova feedback efficiency (i.e., e = 0.25). Fig. 11 shows that for e = 0.5 and 1, for the 
SN rate value oi rj/rjc = 0.1, the velocity dispersions of the H I gas are of the order ~ 5 and 
~ 6 km s -1 , respectively, in closer agreement with the observations. 

On the other hand, the velocity dispersions observed in NGC 628 and NGC 6946 (same 
for the observations presented in Fig. 2), in addition of being affected to some degree by beam 
smearing due to the limited resolution of radio telescopes (this is a minor effect at the outer 
galactic radii), are global velocity dispersions which do not disentangle the true dynamical 
value, which is the only one we have in our simulations, from the thermal broadening of 
the line, simply because the temperature structure of the gas in the observations is not 
known. Thermal broadening, vt = (2 h T/m) 1 ^ 2 , is the dispersion of the velocity probability 
distribution function for an ensemble of particles of a non relativistic and non degenerate 
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gas which is in thermal equilibrium at a kinetic temperature T, where m is the mass of the 
particles and in our case, m is the equal to the proton mass. Seeking a better match between 
our models and the observations, we have corrected for the effect of thermal broadening in 
two ways : 

(a) the simplistic method : We can subtract, quadratically, from the observational 
velocity dispersion values, a velocity component associated with thermal broadening at a 
given equilibrium temperature for the Hi gas, (Jdyn = (<7 2 ot — ff,) 1 / 2 , where a^y-n is the 
dynamical component that is compared with the dynamical component observed in the 
simulations and a b. s and vt are the observed velocity dispersion and the above mentioned 
thermal component, respectively. This is obviously a simplification as we assume the Hi gas 
to have the same temperature at all radii. Fig. 19 shows the corrected velocity dispersions, 
using this method, for three values of the equilibrium temperature of the Hi gas, namely 
100 K, 500 K, and 2000 K. The correction improves the agreement between the observations 
and the simulations, particularly at the outer galactic radii and if the Hi is assumed to be 
warmer than 100 K (bottom plots in Fig. 19 where the Hi gas in NGC 628 and NGC 6946 is 
assumed to be at 2000 K). The fact that the gas in the outer parts could have temperatures 
higher than 100 K is plausible since the density in the outer parts is low and the gas could 
be more easily heated by cosmic rays and photoelectric effect. This is unlikely however 
to be the case for the inner parts of the galaxy where the Hi gas is more likely to have, 
like in the Milky Way, a non-negligible fraction of the H I in the cold phase at around 100 
K. Nevertheless, recent observations by Heiles (2001) suggest that about half of the mass 
of the diffuse interstellar gas in the Galaxy may have temperatures which are larger than 
100 K (a few hundreds to a few thousands Kelvin). Only an accurate determination of 
the temperature structure in galaxies such as NGC 628 and NGC 6946 may lay out strong 
constraints on the contribution of thermal broadening as a function of radius to the total 
velocity dispersion. 

(b) the less simplistic method : Here, we have assumed that the particles in each 
cell have a Gaussian velocity profile which is centered around the local dynamical velocity 
and which have a dispersion in the velocity space equal to vt, where vt is the local thermal 
broadening calculated using the local temperature. The amplitude of the profile is given by 
the local density. Only cells with T < 12000 K and n > 0.25 cm -3 are taken into account. 
The individual velocity profiles are summed up in the velocity space and binned with a 
spectral bin size of 1 km s -1 . Fig. 20 and Fig. 21 display two mass weighted, thermally 
broadened velocity profiles corresponding to simulations (rj/rja, e)=(l, 0.25) and (0.1,0.25), 
respectively. They can be compared to Fig. 8 and Fig. 10, respectively, in order to appreciate 
the effects of thermal broadening on the line profile. Fig. 22 shows a comparison of the 
Hi velocity dispersion ohi to the observational data using this more reliable approach for 
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correcting for the effect of thermal broadening. The result is somewhat encouraging. For 
SFR/Area in the range of 5 x 1CT 3 — 1CT 2 M yr -1 kpc -2 (i.e., rj/rja — 0.5 — 1), the agreement 
to the observations concerning both NGC 628 and NGC 6946 is quite acceptable. At the 
lower SN rates values, the agreement is less satisfying even when the Kennicutt rate-adapted 
average densities are used (open square and star in Fig. 22). However, a comparison only to 
the data of NGC 628 and NGC 6946 might be slightly misleading. In the case of galaxies 
such as NGC 1058 (Dickey, Hanson & Helou 1990) and NGC 3938 (van der Kruit & Shostak 
1982), for which we unfortunately do not have radially dependent estimates of the star 
formation rate, the velocity dispersion levels off in the outer radii at a value of the order of 
5-6 km s" 1 which is in better agreement to the values coming out from our simulations. 



8. THE NEED FOR IMPROVED MODELS 

Effect of the background heating : An important effect, which is not accounted for 
in our models is a background heating of the gas by the photoelectric effect, cosmic rays and 
soft X-rays (see also §2). These background heating processes might be of little importance 
in the case of media with SFR/Area > 5 x 10 -3 M yr _1 kpc -2 (i.e., rj/rjo > 0.5), but 
might play a significant role in maintaining a warmer phase of the gas at the lower SN 
rates. Thus, the velocity line profiles might be broader than what we have calculated in the 
absence of such processes. Kim (2004) showed that the velocity dispersion of the gas in a 
SN driven medium is reduced if the average pressure of the gas in increased. In the absence 
of background heating, the gas cools efficiently to the minimum temperature of 100 K in 
the regions where SN explosions are rare (i.e., this is particularly true for the simulations 
with the low SN rates). Hence, the gas pressure is reduced and the velocity dispersion 
enhanced. On the other hand, the shock produced by a SN explosion expanding in a cold 
medium would lead to higher compressions than a shock propagating in a warmer medium. 
The stronger compressions would lead to an enhanced cooling in the compressed SN shells 
which causes the energetic content of the explosion to be depleted faster. The intensity of 
the background heating is difficult to estimate for systems with different SN rates. As no 
established formulations of this problem exist in the literature, we intend, in future work, to 
model the background heating as being a fraction of the total SN heating and quantify its 
effects on the resulting velocity dispersion of the gas. 

Effect of the vertical structure : The vertical stratification, which we have neglected 
in this work, might be one of the physical effects that we need to include first in subsequent 
models. We intend to perform models with a much larger length of the box in the vertical 
direction, using outflow boundary conditions to allow the gas to escape from the upper and 
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lower boundaries of the simulation box. The escaping hot gas would not affect the Hi 21 
cm line profile. However, denser blobs of gas which are expelled at higher latitudes by 
the SN explosions will cool and fall back into the disk in the form of high velocity clouds. 
In principle, these high velocity clouds should be observed as broad wings in the velocity 
profile and are usually fitted with a second, broader Gaussian function. However, in-falling 
decelerated gas, close to the galactic disk will be most probably mixed with the local gas 
in the disk, thereby contributing, eventually, to the broadening of the velocity profile. The 
impact of HVCs on the galactic can also substantially enhance the local level of turbulence 
particularly for massive HVCs. 

Effect of the chemistry : In the simulations presented in this work, the cooling 
curve we have used is a solar metallicity curve that assumes chemical equilibrium. This 
cooling curve describes principally the radiative cooling by neutral atoms whereas molecular 
cooling by molecules such as the H 2 molecule is neglected. The additional cooling at lower 
temperatures might enhance the local effects of TI and increase the value of the velocity 
dispersion. At the high and intermediate SN rates, this might not be of much relevance as 
the over-densities that are produced are a factor 3-5 the average density (a few cm" 3 ) (see 
Fig. 5) and no dense molecular material is expected to form. At low rates, clouds have time to 
form and condense further before being destroyed by the next generation of SNe (see Fig. 6). 
Molecular hydrogen, starts to form when densities of the order of 10 3 cm" 3 are reached and 
which become shielded against UV radiation (Bergin et al. 2004). Such densities are not 
reached in our simulations, essentially because of the limitations in the numerical resolution. 
The introduction of a simple chemical network to follow locally the fraction of molecular 
hydrogen in a hydrodynamical adaptative-mesh-refinement (AMR) code would help tackle 
the problem of molecules formation in the expanding shells more accurately, and help better 
account for the additional molecular cooling and its effects on the velocity field. 

Effect of the metallicity : The occurrence and efficiency of the thermal instability is 
intimately related to the shape of the cooling curve which characterizes the medium. Cooling 
curves are a reflection of the strength of the emission lines of atoms present in the medium. 
At lower metallicity, the emission lines from metals are weaker and the cooling less efficient. 
The dependence of the cooling rate on the metallicity has been calculated by Boehringer & 
Hensler (1989), unfortunately only in the temperature range 10 4 — 10 8 K. Their results show 
that the cooling rates may differ, in some temperature ranges, by several orders of magnitude 
for metallicities between 10~ 2 < Z/Z Q < 2 (Z/Z Q = 2 is the upper metallicity limit in their 
calculations) and becomes independent of the metallicity for Z/Z Q < 10~ 2 . In future work, 
we plan to investigate the role of metallicity on the dynamics of the ISM by including it as 
a parameter in our simulations. 
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Effect of self-gravity : Kim et al. (2003) in some of their MRI simulations which 
include gravity show that for a single-phase medium and for a Toomre parameter Q th > 1.7 
which is more appropriate for the external regions of galactic disks where the gas surface 
density is low, self-gravity can be responsible for only 20 % of the level of 1.6 — 3.2 km s _1 
turbulence generated by the MRI. However, in our models with SN explosions, the presence 
of self-gravity, provided enough numerical resolution is affordable to resolve the expanding 
shells, might lead to the development of Rayleigh- Taylor and Kelvin-Helmotz instabilities 
(Shu 1992), particularly in the case of SN remnants evolving in non-spherical environments 
(see e.g., Gazol-Patino & Passot 1999). 

Effect of the magnetic field : In the present work, magnetic fields have been ne- 
glected. Several authors have performed simulations of a supernova explosions in a magne- 
tized medium (e.g., Ferrierre 1991; Gazol-Patino & Passot 1999; Korpi et al. 1999a,b; Kim 
2004; de Avillez & Breitschwerdt 2005; Mac Low et al. 2005). The latter simulations agree 
that the effect of magnetic fields in essentially to oppose the radial expansion of a supernova 
remnant, thus reducing the energy transmitted to the ISM and reducing the velocity disper- 
sion. In particular, the results of Kim (2004) show that the velocity dispersion in a SN driven 
medium is related to the total (thermal + magnetic) pressure. Dispersions are smaller in en- 
vironments with higher total pressures (i.e., higher magnetic field values for a given thermal 
pressure). For a comparison, at the Galactic SN rate we find that the one- dimensional total 
gas velocity dispersion is ~ 8.5 km s _1 , whereas in the magnetized case Kim (2004) finds ~ 9 
and ~ 10 km s" 1 for a weak magnetic field value of 2 /xG associated to an average density of 
0.2 cm -3 in the directions parallel to the mean field and perpendicular to it, respectively. For 
a stronger field value of 8 /zG associated to an average density of 0.8 cm~ 3 , Kim (2004) finds 
velocity dispersions of ~ 5 and ~ 7 km s _1 in the directions parallel and perpendicular to 
the field, respectively. In view of Kim's (2004) results, we speculate that in our simulations, 
as SNe occur in region of higher density, the existence of a magnetic field which would be 
compressed in those regions, would have the effect of lowering our measured values of the 
total gas velocity dispersion by a factor of a few tens of percent (typically 10-20 %). 



9. SUMMARY AND DISCUSSION 

In this paper, we investigated the dependence of the velocity dispersion in the interstellar 
medium (ISM) on the supernova (SN) rate 77, the SN feedback efficiency e, and in some cases 
on the ISM average number density n. We use local, three-dimensional numerical simulations 
in which SN type II explosions are detonated in random positions of the grid, separated 
by time intervals which are inversely proportional to the SN rate. Radiative cooling of 
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the gas is also taken into account with a minimum cutoff temperature of 100 K. For the 
purpose of simplifying the problem, other physical processes and characteristics of galactic 
disks such as the vertical stratification, magnetic fields and gravity are neglected. For each 
simulation, we calculate the three-dimensional characteristic velocity dispersion v c (Eq. 1) 
and the one-dimensional velocity dispersion, a, obtained by fitting the line of sight velocity 
profile with a Gaussian function. We also calculate am, which is the one-dimensional velocity 
dispersion obtained from a Gaussian fit of the line of sight velocity profile of the gas with a 
temperature T < 12000 K and a number density n > 0.25 cm -3 (i.e., the Hi gas). Ideally, a 
full investigation of the two-dimensional parameter space (77, e) would be necessary, however, 
this would lead us beyond our current computational capabilities. At this stage, we have 
taken an intermediate approach and explored the effects of 77 and e independently by fixing 
one parameter and varying the other. 

Our results show that v c , a and am depend on the SN feedback efficiency e as Ay/e, 
where A is a scaling coefficient, different for each quantity. This is expected as the velocity 
has the dimensions of the square root of an energy. In a second set of simulations, we 
fixed e to a value of 0.25 and varied the SN rate in the range [0.01,10] r)a, where r)o is the 
Galactic SN rate. We compare the velocity dispersion for the simulations with the different 
SN rates to NGC 628 and NGC 6964 for which both the velocity dispersion and the the 
star formation rates are radially resolved. The dependence of the velocity dispersion on the 
SN rate is complex. For values of r]/r] G > 0.5 (i.e., SFR/Area > 10~ 2 M yr^ 1 kpc~ 2 ), v c 
and a increase sharply with increasing values of rj. For rj/r/c < 0.5, v c and a show a slower 
decrease with a decreasing 77. This transition at r]/r] G ~ 0.5 is probably an indication that 
the velocity dispersion of the medium in the low rate regime is not fixed by SN driving alone. 
We interpret the flatness of the a — rj relation as resulting from the efficient development of 
thermal instability (TI) in the medium with the low SN rates. We quantify the efficiency of 
TI by evaluating the volume filling factor of the unstable gas F u for the different SN rates. 
F u appears to be correlated with the velocity dispersion. For large values of the SN rate 
Tj/rjc > 0.5 — 1, F u is close to zero, with a transition to non zero values at rj/r/c ~ 0.5, peaks 
at rj/r/c ~ 0.1, before decreasing slowly at smaller SN rates. Interestingly, the position of the 
transition to the starburst regime (high SFR rates associated with high velocity dispersions) 
is in relatively good agreement with the one seen in the observations (see Fig. 2, Fig. 18 and 
Fig. 19). 

We compare am to the observations as a function of the SN rate (i.e., SFR rates), a 
and am have nearly similar values at the low SN rates as most of the gas in this regime is 
Hi gas. Hi gas is not present in the simulations with ri/ric > 1- am is observed to be nearly 
independent of the SN rate, leveling off at ~ 3 km -1 . This value is a factor 2 — 3 lower that 
the velocity dispersion plateau observed in galaxies such as NGC 1058, NGC 628 and NGC 
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6964 (Fig. 1 and Fig. 18). However, the set of simulations with r]/i] G = 0.1 and the different 
feedback efficiencies (i.e, Fig. 11) suggests that Hi gas velocity dispersions of the order of 
5-6 km s" 1 can be obtained for feedback efficiencies e > 0.5, in closer agreement with the 
observations. For the simulations with our fiducial value of e, we have corrected the H I line 
velocity profiles by accounting for the effect of thermal broadening. A reasonable agreement 
is found between the models and the observations for values of the SFR/Area in the range 
5 x 10~ 3 - 10~ 2 M yr _1 kpc~ 2 . For smaller values of the SFR, the fact that the velocity 
dispersions are a factor ~ 2 smaller than the observed values could result from the fact that 
we have underestimated the SN feedback efficiency. Otherwise, it might be an indication of 
the existence of secondary heating and/or driving mechanisms in the outer parts of galaxies 
where the star formation rate is low. An investigation of the effects of SN driving in the 
presence of other physical processes and gas instabilities is essential and is left to future 
work. 

It is a pleasure to thank Enrique Vazquez-Semadeni, Fabian Walter, Hans- Walter Rix, 
Axel Brandenburg, Maarit Korpi, Joop Schaye and Robert Piontek for useful comments and 
discussion. We are also very grateful to the anonymous referee for interesting comments 
and suggestions. S. D. would like to thank Hans- Walter Rix and Thomas Henning for the 
financial support at the MPIA. Calculations have been performed on the MPIA's SGI Origin 
2000 located at the Rechnenzentrum of the Max-Planck Gesellschaft, Garching. ZEUS-3D 
was used by courtesy of the Laboratory of Computational Astrophysics at the NCSA. This 
research has made use of NASA's Astrophysics Data System Bibliographic Services. 

REFERENCES 

Abbott, D. C. 1982, ApJ, 263, 723 

Avila- Reese, V., & Vazquez-Semadeni, E. 2001, ApJ, 553, 645 
Balbus, S. A., & Hawley, J. F. 1991, ApJ, 376, 214 

Bergin, E. A., Hartmann, L. W., Raymond, J. C, & Ballesteros-Paredes, J. 2004, ApJ, 612, 
921 

Boehringer, H., & Hensler, G. 1989, A&A, 215, 147 
Boroson, T. 1981, ApJS, 46, 177 



-19- 



Bureau, M., Walter, F., van Gorkom, J., & Carignan, C. 2004, in Recycling intergalactic 
and interstellar matter, IAU Symposium series, Vol. 217, ed. P.-A. Due, J. Braine, & 
E. Brinks 

Cappellaro, E. Evans, R., & Turatto, M. 1999, A&A, 351, 459 

Chevalier R. A. 1977, ARA&A, 15, 175 

Dalgarno, A., & McCray, R. A. 1972, ARA&A, 10, 375E 

de Avillez, M. A. 2000, MNRAS, 315, 479 

de Avillez, M. A., & Berry, D. 2001, MNRAS, 328, 708 

de Avillez, M. A., & Breitschwerdt, D. 2004, A&A, 425, 899 

de Avillez, M. A., & Breitschwerdt, D. 2005, A&A, 436, 585 

de Blok, W. J. G., & Walter, F. 2005, AJ, accepted,(astro-ph/0509107) 

de Jager, G., & Davies, R. D. 1971, MNRAS, 153, 9 

Dib, S., & Burkert, A. 2005, ApJ, 630, 238 

Dickey, & J. M. Lockman, F. J. 1990, ARA&A, 28, 215 

Dickey, J. M., Murray-Hanson, M., & Helou, G. 1990, ApJ, 352, 522 

Dopita, M. A., & Ryder, S. D. 1994, ApJ, 430, 163 

Dziourkevitch, N., Elstner, D., & Riidiger G., 2004, A&A, 423, 29 

Dziourkevitch, N. 2005, in The magnetized Plasma in Galaxy Evolution, ed. K. T. Chyzy, K. 
Otmianoswka-Mazur, M. Soida, & R.-J. Dettmar (Jagiellonian University, Krakow), 
p 74 

Elmegreen, B., & Scalo, J. 2004, ARA&A, 42, 211 
Elmegreen, D. M., Elmegreen, B. G. 1984, ApJS, 54, 127 
Ferriere, K. 1992a, ApJ, 389, 286 
Ferriere, K. 1992b, ApJ, 391, 188 
Ferriere, K. 1998a, ApJ, 497, 759 



-20- 

Ferriere, K. 1998b, ApJ, 503, 700 

Ferriere, K., Mac Low, M.-M., & Zweibel, E. G. 1991, ApJ, 375, 239 
Ferriere, K., & Schmitt, D. 2000, A&A, 358, 125 

Ferguson, A. M. N., Wyse, R. F. G., Gallagher, J. S., & Hunter, D. 1998, ApJ, 506, 19 

Fioc, M., & Rocca-Volmerange, B. 1997, A&A, 326, 950 

Gazol-Patino, A., & Passot, T. 1999, ApJ, 519, 748 

Gordon, K. J. 1971, ApJ, 169, 235 

Heiles, C. 1987, ApJ, 315, 555 

Heiles, C. 2001, ApJ, 551, L105 

Hippelein, H., Haas, M., Tuffs, R. J., Lemke, D., Stickel, M., Klaas, U., & Volk, K. J. 2003, 
407, 137 

Huchtmeier, W. K. 1973, A&A, 22, 91 

Kamphuis, J., & Sancisi, R. 1993, A&A, 273, L31 

Kennicutt, R. C. 1998, ApJ, 498, 541 

Kennicutt, R. C. 1998, ARA&A, 36, 189 

Kessel-Deynet, O., & Burkert, A. 2003, MNRAS, 338, 545 

Kim, J., Balsara, D., & Mac-Low, M.-M. 2001, JKAS, 34, 333 

Kim, J. 2004, JKAS, 37, 237 

Kim, W.-T., Ostriker, E. C, & Stone, J. M. 2003, ApJ, 599, 1157 

Korpi, M. J., Brandenburg, A., Shukurov, A., Tuominen, L, & Nordlund, A. 1999, ApJ, 
514, 9 

Korpi, M. J., Brandenburg, A., Shukurov, A., & Tuominen, I. 1999, A&A, 350, 230 
Larson R. B. 1981, MNRAS, 194, 809 

Li, Y., Mac Low, M.-M., & Klessen, R. S. 2005, ApJ, 626, 823 

Li, Y., Mac Low, M.-M., & Klessen, R. S. 2005, ApJ, submitted,(astro-ph/0508054) 



- 21 - 



Martin, C. L., & Kennicutt, R. C. 2001, ApJ, 555, 301 

Mac Low, M.-M., & Klessen, R. 2004, Rev. Mod. Phys., 76, 125 

Mac Low, M.-M., Balsara, D. S., Kim, J., & de Avillez, M. A. 2005, ApJ, 626, 864 

Meurer, G. R., Mackie, G., & Carignan, C. 1994, AJ, 107, 2021 

Meurer, G. R., Carignan, C, Beaulieu, S. F., & Freeman, K. C. 1996, AJ, 111, 1551 

Miller, G. E., & Scalo, J. M. 1979, ApJS, 41, 513 

Mckee, C. F., & Ostriker, J. P. 1977, ApJ, 218, 148 

Piontek, R. A., & Ostriker, E. C. 2004, ApJ, 601, 905 

Piontek, R. A., & Ostriker, E. C. 2005, ApJ, 629, 849 

Prantzos, N., & Silk, J. 1998, ApJ, 507, 229 

Riccoti, M., & Ferrara, A. 2002, A&A, 334, 684 

Rohlfs, K., & Wilson, T. L. 1996, Tools of Radio Astronomy (Berlin: Springer) 

Rosen, A., & Bregman, J. N. 1995, ApJ, 440, 634 

Rownd, B. K., Dickey, J. M., & Helou, G. 1994, ApJ, 108, 1683 

Salpeter, E. E. 1955, ApJ, 121, 161 

Scalo, J. M. 1987 in Interstellar Processes, ed. D. J. Hollenbach & H. A. Thronson (Dor- 
drecht :Reidel), 349 

Schartmann, M., Meisenheimer, K., Camenzind, M., Wolf, S., & Henning, T. 2005, A&A, 
437, 861 

Schaye, J. 2004, ApJ, 609, 667 

Schmidt, M. 1959, ApJ, 129, 243 

Schmidt, M. 1963, ApJ, 137, 758 

Selwood, J. A., & Balbus, S. A. 1999, ApJ, 551, 660 

Shostak, G. S., & van der Kruit, P. C. 1984, A&A, 132, 20 

Shu, F. H. 1992, Physics of Astrophysics, Vol. II (Mill Valley : University Science Books ) 



-22 - 



Slyz, A. D., Devriendt, J. E. G., Bryan, G., & Silk, J. 2005, MNRAS, 356, 737 

Stil, J. M., & Israel, F. P. 2002, A&A, 392, 473 

Stone, J. M., & Norman, M. L. 1992a, ApJS, 80, 753 

Stone, J. M., & Norman, M. L. 1992b, ApJS, 80, 791 

Sutherland, R. S., & Dopita, M. A. 1993, ApJS, 88, 253 

Tammann, G. A., Iroffler, W., & Schroder, A. 1994, ApJS, 92, 487 

Thornton, K., Gaudlitz, M., Janka, H.-Th., & Steinmetz, M. 1998, ApJ, 500, 119 

van den Bergh, S., & McClure, R. D. 1990, ApJ, 359, 277 

van der Kruit, P. C, & Shostak, G. S. 1982, A&A, 105, 351 

van Zee, L., Haynes, M., Salzer, J. J., & Broeils, A. H. 1997, AJ, 113, 5 

van Zee, L. 2001, AJ, 121, 2003 

Vazquez-Semadeni, E., Gazol, A., & Scalo, J. 2000, ApJ, 540, 271 
Wada, K., & Norman, C. A. 1999, ApJ, 516 , 13-16 
Wada, K., & Norman, C. A. 2001, ApJ, 547, 172 
Wada, K., Spaans, M., & Kim, S. 2000, ApJ, 540, 797 
Wada, K., Meurer, G., & Norman, C. A. 2002, ApJ, 577, 197 
Woosley, S. E., & Weaver, T. A. 1986, ARA&A, 24, 205 



This preprint was prepared with the AAS IATgX macros v5.2. 



-23- 




Fig. 1. — Radial dependence of the velocity dispersion in NGC 1058. Adapted from Dickey, 
Hanson & Helou (1990) (Fig. 3 in their paper), assuming a distance of 10.2 Mpc (Boroson 
1981). 
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Fig. 2. — Characteristic velocity dispersion a of a sample of galaxies as a function of the 
surface averaged star formation rate (units of M yr -1 kpc~ 2 ). The latter value is calculated 
from H a observations. The galactic surface is calculated as being n (3 r^) 2 , where is the 
radial length scale of each galaxy. References are (1) van Zee et al. (1997) (2) van Zee (2001) 
(3) Meurer et al. (1996) (4) Hutchmeier (1973) (5) Gordon (1971) (6) de Jager & Davies 
(1971) (7) Hippelein et al. (2003) (8) Elmegreen & Elmegreen (1984) (9) Meurer, Mackie & 
Carignan (1994) (10) Stil & Israel (2002). 
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Fig. 3. — Convergence study of the shell radius for a supernova explosion occurring in 
medium with initial temperature of 10 4 K and uniform and homogeneous density of 0.5 cm -3 . 
The size of the box is 1 kpc. The three resolutions correspond to 64 3 (green, hexagons), 128 3 
(blue, squares), and 256 3 (red, triangles). 
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Fig. 4. — Time evolution of the kinetic energy in a number of selected simulations. The 
supernova rate, normalized to the Galactic value is shown for each curve. The feedback 
efficiency for all three models is e = 0.25. 
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Fig. 5. — Density cut in the mid-plane of the data cube at time=72 Myrs for the model with 
(r)/r)G, e) = (0.1,0.25). The color bar data scale is logarithmic. 
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Fig. 6. — Density cut in the mid-plane of the data cube at time=60 Myrs for the model with 
{v/Vg, e ) — (1; 0.25). The color bar data scale is logarithmic. 
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Fig. 7. — Mass weighted velocity profile for the model with (t)/t]q, e) = (l,0.25) and a spectral 
bin size of 1 km s _1 . The profile samples all gas elements in the simulation box. Over-plotted 
is the Gaussian fit (dotted line). 
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Fig. 8. — Mass weighted velocity profile for the model with (rj/rj G , e) = (l,0.25) and a spectral 
bin size of 1 km s _1 . The profile samples only gas elements with T < 12000 K and n > 0.25 
cm -3 . Over-plotted is the Gaussian fit (dotted line). 
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Fig. 9. — Mass weighted velocity profile for the model with (rj/rja, e)=(0. 1,0.25) and a 
spectral bin size of 1 km s _1 . The profile samples all gas elements in the simulation box. 
Over-plotted is the Gaussian fit (dotted line). 
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Fig. 10. — Mass weighted velocity profile for the model with (77/770, e)=(0.1,0.25) and a 
spectral bin size of 1 km s -3 . The profile samples only gas elements with T < 12000 K and 
n > 0.25 cm -3 . Over-plotted is the Gaussian fit (dotted line). 
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Fig. 11. — Characteristic velocity v c (Eq. 1) (triangles), mass- weighted line of sight velocity a 
(diamonds) and H I gas mass- weighted line of sight velocity cthi (squares) (both have spectral 
bin size of 0.1 km s~ x ) as a function of the supernova feedback energy equal to e x 10 51 ergs. 
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Fig. 12. — Characteristic velocity v c as a function of supernova rate 77 (normalized to the 
Galactic value i] G ). 
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Fig. 13. — Velocity dispersion of the gas a as a function of the supernova rate r\ (normalized 
to the Galactic value t]g, filled triangles). Filled square and filled hexagon correspond to 
simulations where the average density has been reduced by factors of 5 and 10, respectively. 
The displayed data has been obtained using a spectral bin size of 1 km s -1 . The constant 
line at 6 km s -1 is drawn to guide the eye. 
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Fig. 14. — Volume filling factor of the unstable gas (398 K < T < 10000 K) as a function 
of the supernova rate 77 (normalized to the Galactic value rja)- The error bars represent the 
statistical errors for 5 estimates of each value (last 5 Myrs in each simulation). 
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Fig. 15. — Time evolution of the volume filling factor of the unstable gas (398 K < T < 1000 
K, dashed line), cold gas (T < 398 K, solid line) and warm gas (T > 10000 K, dot-dashed 
line) for the simulation with the set of parameters (rj/rja, e) = (0.01, 0.25). 



- 38 - 



u 

- 0.4 

0.2 
0.0 



1 , 
















p > ' 

T\ i 








t 1 '' 


\ 

\ 

\ 

\ 

\ 

\ 

/ ^ - ~~ ~ ~ - - 






J- i; 

+ ( 
I 

f M 






J ' ' 








/ 1 








1 1 








\ 












I 





20 40 00 80 

time (in Myrs) 



Fig. 16. — Time evolution of the volume filling factor of the unstable gas (398 K < T < 1000 
K, dashed line), cold gas (T < 398 K, solid line) and warm gas (T > 10000 K, dot-dashed 
line) in the simulation with the set of parameters (r]/r]G,€) = (0.1,0.25). 
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Fig. 17. — Velocity dispersion of the gas a as a function of the star formation rate per unit 
area (M yr -1 kpc~ 2 ) (filled squares), open square and star correspond to simulations where 
the density has been reduced by a factor of 5 and 10, respectively. The displayed data has 
been obtained using a spectral bin size of 1 km s -1 . For NGC 628 (open diamonds) and NGC 
6946 (open triangles), the velocity dispersion is derived from Hi 21 cm line observations and 
the star formation rates from Ha observations. 
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Fig. 18. — Velocity dispersion of the Hi gas er H i as a function of the star formation rate per 
unit area (M yr -1 kpc~ 2 ) (filled squares), open square and star correspond to simulations 
where the density has been reduced by a factor of 5 and 10, respectively. The displayed data 
has been obtained using a spectral bin size of 1 km s" 1 . For NGC 628 (open diamonds) and 
NGC 6946 (open triangles), the velocity dispersion is derived from H I 21 cm line observations 
and the star formation rates from Ha observations. 
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Fig. 19. — Same as in Fig. 18 but where thermal broadening has been subtracted from the 
observed velocity dispersion assuming that the Hi gas is at 100 K (top), 500 K (middle) and 
2000 K (bottom). Spectral bin size used is 1 km s _1 . 
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Fig. 20. — Mass weighted velocity profile for the model with (rj/rjc, e) = (l,0.25) and corrected 
for thermal broadening. The spectral bin size of 1 km s _1 . The profile samples only gas 
elements with T < 12000 K and n > 0.25 cm' 3 . Over-plotted is a Gaussian fit (dotted line). 
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Fig. 21. — Mass weighted velocity profile for the model with (rj/rja, e)=(0. 1,0. 25) and cor- 
rected for thermal broadening. The spectral bin size of 1 km s _1 . The profile samples only 
gas elements with T < 12000 K and n > 0.25 cm -3 . Over-plotted is a Gaussian fit (dotted 
line) . 
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Fig. 22. — Velocity dispersion of the Hi gas a H , asa function of the star formation rate per 
unit area (M© yr" 1 kpc" 2 ) (filled squares). The values have been corrected for the effect of 
thermal broadening. Open square and star correspond to simulations where the density has 
been reduced by a factor of 5 and 10, respectively. The displayed data has been obtained 
using a spectral bin size of 1 km s -1 . 



